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Abstract.  This  paper  studies  Lagrangian  mixing  in  a  two- 
dimensional  barotropic  model  for  hurricane-like  vortices. 
Since  such  flows  show  high  shearing  in  the  radial  direction, 
particle  separation  across  shear-lines  is  diagnosed  through  a 
Lagrangian  field,  referred  to  as  R-field,  that  measures  trajec¬ 
tory  separation  orthogonal  to  the  Lagrangian  velocity.  The 
shear-lines  are  identified  with  the  level-contours  of  another 
Lagrangian  field,  referred  to  as  S-field,  that  measures  the 
average  shear-strength  along  a  trajectory.  Other  fields  used 
for  model  diagnostics  are  the  Lagrangian  field  of  finite-time 
Lyapunov  exponents  (FTLE- field),  the  Eulerian  Q- field,  and 
the  angular  velocity  field.  Because  of  the  high  shearing,  the 
FTLE- field  is  not  a  suitable  indicator  for  advective  mixing, 
and  in  particular  does  not  exhibit  ridges  marking  the  location 
of  finite-time  stable  and  unstable  manifolds.  The  FTLE- field 
is  similar  in  structure  to  the  radial  derivative  of  the  angu¬ 
lar  velocity.  In  contrast,  persisting  ridges  and  valleys  can 
be  clearly  recognized  in  the  R-field,  and  their  propagation 
speed  indicates  that  transport  across  shear-lines  is  caused  by 
Rossby  waves.  A  radial  mixing  rate  derived  from  the  R- 
field  gives  a  time-dependent  measure  of  flux  across  the  shear¬ 
lines.  On  the  other  hand,  a  measured  mixing  rate  across  the 
shear-lines,  which  counts  trajectory  crossings,  confirms  the 
results  from  the  R -field  mixing  rate,  and  shows  high  mixing 
in  the  eyewall  region  after  the  formation  of  a  polygonal  eye- 
wall,  which  continues  until  the  vortex  breaks  down.  The  lo¬ 
cation  of  the  R -field  ridges  elucidates  the  role  of  radial  mix¬ 
ing  for  the  interaction  and  breakdown  of  the  mesovortices 
shown  by  the  model. 
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1  Introduction 

Several  recent  studies  (Frank  and  Ritchie,  1999,  2001;  Mont¬ 
gomery  et  al.,  2006;  Hendricks  and  Schubert,  2009)  are  de¬ 
voted  to  the  mixing  of  fluid  from  different  regions  of  a  hurri¬ 
cane,  which  is  considered  as  a  fundamental  mechanism  that 
is  intimately  related  to  hurricane  intensity.  A  complete  un¬ 
derstanding  of  these  mixing  processes,  in  particular  the  eye- 
eyewall  mixing  (Bryan  and  Rotunno  2009a,  b;  Cram  et  al., 
2007;  Braun  et  al.,  2006;  Montgomery  et  al.,  2002,  2006; 
Willoughby,  2001),  is  expected  to  improve  our  understanding 
of  the  physical  mechanisms  that  regulate  hurricane  intensity. 
Since  mixing  is  based  on  particle  motion,  the  Lagrangian 
frame  of  reference  provides  the  most  natural  framework  in 
which  it  can  be  diagnosed.  Much  progress  has  been  made 
in  recent  years  in  the  study  of  Lagrangian  mixing  in  two- 
dimensional  incompressible  flows  (Haller  and  Poje,  1997; 
Haller  and  Yuan,  2000;  Haller,  2000,  2001,  2002;  Shadden 
et  al.,  2005),  resulting  in  a  number  of  different,  though  re¬ 
lated  diagnostics,  most  of  which  are  based  on  concepts  from 
dynamical  systems  theory.  For  applications  of  Lagrangian 
techniques  to  atmospheric  models,  see  Joseph  and  Legras 
(2001)  and  Huber  et  al.  (2001). 

Much  insight  into  specific  aspects  of  mixing  in  hurricanes 
can  be  gained  from  the  study  of  simplified  two-dimensional 
models.  Basically  there  are  two  classes  of  such  models:  Ax- 
isymmetric  models,  most  notably  the  model  of  Rotunno  and 
Emanuel  (1987),  and  planar  models  such  as  the  model  of 
Kossin  and  Schubert  (2001)  and  Schubert  et  al.  (1999).  In 
this  paper  we  apply  Lagrangian  techniques  to  analyze  mix¬ 
ing  in  the  planar,  nondivergent  barotropic  model  of  Kossin 
and  Schubert  (2001).  Our  analysis  confirms  a  study  of 
Kossin  and  Eastin  (2001)  which  illustrates  that  significant 
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eye-eyewall  mixing  occurs  during  polygonal  eyewall  tran¬ 
sitions.  Lagrangian  mixing  in  the  axisymmetric  model  in¬ 
troduced  in  Rotunno  and  Emanuel  (1987)  is  investigated  in 
Rutherford  et  al.  (2009).  For  another  discussion  of  axisym¬ 
metric  mixing,  see  Wirth  and  Dunkerton  (2006). 

The  model  studied  in  this  paper  provides  a  two- 
dimensional  representation  of  a  hurricane  that  initiates  with 
an  annular  ring  of  enhanced  vorticity,  and  then  undergoes  a 
vortex  breakdown  resulting  in  a  monopolar  end  state.  Dur¬ 
ing  the  breakdown,  a  polygonal  eyewall  occurs,  which  forms 
four  elliptical  pools  of  high  vorticity.  Mixing  of  potential 
vorticity,  which  in  this  model  is  proportional  to  relative  vor¬ 
ticity,  can  be  visualized  using  Eulerian  diagnostic  measures 
of  instantaneous  particle  separation.  A  commonly  used  Eu¬ 
lerian  diagnostic  is  the  so  called  Q -field,  derived  from  the 
Jacobian  of  the  Eulerian  velocity  field.  According  to  the 
Okubo-Weiss  criterion  (Schubert  et  al.,  1999),  positive  val¬ 
ues  of  this  field  indicate  instantaneous  particle  separation, 
whereas  negative  values  indicate  contraction.  For  our  model, 
the  Q-field  shows  that  regions  of  high  relative  vorticity  gradi¬ 
ent  are  also  places  where  high  trajectory  separation  and  mix¬ 
ing  occurs. 

While  Eulerian  measures  of  mixing  can  only  diagnose  in¬ 
stantaneous  particle  separation,  Lagrangian  techniques  uti¬ 
lize  a  moving  frame  approach  along  trajectories  and  com¬ 
pute  measures  for  the  average  separation  over  a  finite  inte¬ 
gration  time.  This  approach  is  particularly  useful  in  time- 
dependent  velocity  fields,  where  trajectories  may  cross  Eu¬ 
lerian  streamlines  (Dunkerton  et  al.,  2009).  Much  of  the 
recent  work  in  Lagrangian  mixing  has  extended  the  ideas 
of  hyperbolicity  for  steady  flows  to  time  dependent  velocity 
fields  (Haller,  2002;  Haller  and  Poje,  1997;  Haller  and  Yuan, 
2000;  Green  et  al.,  2006;  Ide  et  al.,  2002),  generalizing  the 
concept  of  stable  and  unstable  manifolds  of  an  equilibrium 
to  the  stable  and  unstable  manifolds  of  a  hyperbolic  trajec¬ 
tory.  These  manifolds  are  referred  to  as  Lagrangian  coherent 
structures  (LCS’s).  Even  in  two-dimensional  flows,  time- 
dependence  can  give  rise  to  multiple  intersections  of  these 
manifolds,  leading  to  a  partition  of  the  flow  into  invariant  re¬ 
gions  (lobes),  and  to  mixing  through  the  lobe  dynamics  (Ide 
et  al.,  2002). 

Efficient  visualization  of  LCS’s  is  accomplished  through 
the  construction  of  Lagrangian  scalar  fields,  which  mea¬ 
sure  separation  of  nearby  trajectories.  Current  Lagrangian 
methods  utilize  a  variety  of  fields,  including  finite-time  Lya¬ 
punov  exponents  (Haller,  2002;  Haller  and  Poje,  1997;  Haller 
and  Yuan,  2000),  finite-size  Lyapunov  exponents  (Koh  and 
Legras,  2002;  Green  et  al.,  2006),  and  relative  dispersion 
(Huber  et  al.,  2001).  Each  of  these  methods  defines  a  scalar 
field  and  the  LCS’s  as  maximal  ridges  of  that  field. 

To  study  Lagrangian  mixing  in  our  model,  we  compute 
particle  trajectories  from  the  numerically  calculated,  time- 
varying  velocity  field.  The  Lagrangian  diagnostic  fields  are 
functions  of  the  initial  time  and  position  at  which  the  tra¬ 
jectories  are  seeded.  A  comparison  of  these  fields  with  the 


Okubo-Weiss  criterion  indicates  that  high  particle  separation 
predicted  from  the  Q -field  typically  does  not  coincide  with 
Lagrangian  hyperbolic  structures,  however  the  Lagrangian 
Q -field,  formed  by  integrating  Q -values  along  particle  tra¬ 
jectories,  shows  a  greater  relation  to  other  Lagrangian  fields. 

An  important  feature  of  the  particle  trajectories  calculated 
from  our  model  is  that  they  show  an  almost  circular  motion, 
combined  with  high  shearing  in  the  radial  direction.  The 
problem  caused  by  this  high  shear  is  that  trajectory  separa¬ 
tion  and  mixing  occur  without  the  entrainment  of  trajectories, 
as  the  mixing  is  largely  diffusive.  A  key  question  that  we  aim 
to  answer  is  whether  coherent  structures  that  play  a  role  in  the 
systematic  transport  of  trajectories  can  persist  through  high 
shear. 

Distinct  regions  of  trajectories  with  similar  properties  be¬ 
come  more  difficult  to  distinguish  through  the  use  of  scalar 
fields  which  measure  only  distance,  such  as  the  field  of  finite¬ 
time  Lyapunov  exponents  (TTLL-field).  In  fact,  the  FTLE- 
field  computed  from  our  model  does  not  show  distinguished 
ridges  characteristic  of  hyperbolic  mixing.  Instead,  the  struc¬ 
ture  of  the  TTLL-field  is  very  similar  to  the  structure  of  the 
radial  derivative  of  the  angular  velocity,  indicating  that  the 
TTLL-field  is  dominated  by  the  shear  and  not  by  hyperbolic 
mixing. 

In  order  to  separate  shear  from  hyperbolic  mixing  we  fol¬ 
low  the  approach  used  in  Haller  and  Iacono  (2003),  and  de¬ 
compose  the  separation  of  trajectories  in  the  directions  along 
and  normal  to  the  Lagrangian  velocity.  This  approach  allows 
us  to  identify  two  Lagrangian  fields:  The  R -field,  which  is 
a  diagnostic  for  hyperbolic  mixing  normal  to  the  Lagrangian 
velocity,  and  the  S-field,  which  is  a  measure  of  shearing,  and 
is  used  to  define  shear-lines  by  its  level-contours.  In  contrast 
to  the  TTLE’-field,  the  R -field  shows  distinct  ridges  and  val¬ 
leys  observable  as  coherent  structures.  The  evolution  of  these 
structures  provides  a  mechanism  for  mixing  through  the  eye- 
wall,  and  their  speed  indicates  that  this  mixing  is  caused  by 
Rossby  waves.  The  structures  are  particularly  distinct  after 
polygonal  eyewall  formation,  and  they  persist  until  the  vor¬ 
tex  breaks  down,  in  regions  where  the  Okubo-Weiss  criterion 
predicts  pools  of  high  separation  associated  with  the  forma¬ 
tion  of  pools  of  high  vorticity. 

We  note  that  another  approach  to  diagnosing  mixing  in 
the  presence  of  shear  is  based  on  subtracting  a  mean  shear 
from  the  flow.  This  approach  was  introduced  by  Andrews 
and  McIntyre  (1978)  using  a  generalized  Lagrangian  mean 
for  nonlinear  waves,  and  was  subsequently  developed  further 
and  refined  to  a  modified  Lagrangian  mean  to  quantify  and 
distinguish  stirring  from  irreversible  mixing,  see  McIntyre 
(1980)  and  Dunkerton  (1980). 

The  time-dependence  of  our  velocity  field  leads  to  time- 
dependent  shear-lines,  and  regions  of  high  orthogonal  (hy¬ 
perbolic)  separation  lead  to  sets  of  trajectories  that  are  mixed 
through  the  shear-lines.  We  quantify  this  mixing  by  introduc¬ 
ing  measured  (via  trajectory  counting)  and  predicted  (from 
the  R -field)  mixing  rates.  In  addition,  we  study  radial  mixing 
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rates  defined  by  angular  averages  of  the  TTLE-field,  the  S- 
field,  and  the  R -field.  The  mixing  rates  defined  through  the 
former  two  fields  are  characteristic  of  shearing  and  give  spu¬ 
riously  a  false  sense  of  mixing  during  the  initial  phase  of  the 
model,  where  “true  mixing”  occurs  after  the  polygonal  eye- 
wall  formation. 

In  previous  work  on  the  same  model,  Hendricks  and  Schu¬ 
bert  (2009)  have  applied  the  Lagrangian-Eulerian  hybrid 
method  of  effective  diffusion  (Nakamura,  1996;  Shuckburgh 
and  Haynes,  2003).  Here  diffusive  mixing  properties  are 
computed  based  on  the  increasing  lengths  of  the  vorticity 
contours,  with  the  computations  initialized  at  the  initial  time 
of  the  model.  The  resulting  mixing  rate  is  a  function  of  an 
effective  radius  and  the  integration  time,  and  shows  similar 
structures  as  our  mixing  rates. 

Our  methods  depart  from  those  of  Hendricks  and  Schu¬ 
bert  (2009)  in  that  we  utilize  a  moving  time  window,  which 
attributes  mixing  to  short-time  advective  events.  Rather  than 
determining  contour  lengths,  we  study  transport  across  con¬ 
tours  of  the  S-field.  The  resulting  mixing  rates  are  com¬ 
pletely  determined  by  the  given  velocity  field,  that  is,  they 
do  not  depend  on  a  chosen  initial  profile  of  the  tracer  distri¬ 
bution. 

The  outline  of  the  paper  is  as  follows.  We  begin,  in  Sect.  2, 
with  an  overview  of  the  nondivergent  barotropic  model,  and 
of  the  numerical  methods  used  to  compute  the  velocity  field 
and  the  particle  trajectories.  In  Sect.  3  we  introduce  the  scalar 
fields  utilized  for  diagnosing  mixing  and  shear:  The  Eulerian 
Q -field,  the  Lagrangian  Q -field,  the  angular  velocity  field, 
the  Lagrangian  FTLE- field,  the  R -field,  and  the  S-field.  The 
latter  two  fields  are  extracted  from  the  transformed  varia¬ 
tional  system  introduced  in  Haller  and  Iacono  (2003).  The 
main  results  of  the  paper  are  presented  in  Sects.  4  to  6.  In 
Sect.  4  we  study  the  behavior  of  the  three  Lagrangian  fields 
for  a  fixed  initial  time  of  6  h,  after  a  polygonal  eyewall  has 
formed,  and  for  different  integration  times.  The  ridge,  valley, 
and  edge  structures  observed  in  the  R -field  are  identified  with 
coherent  structures  and  invariant  sets  relative  to  the  shear¬ 
ing.  In  Sect.  5  we  fix  the  integration  time  to  1  h  and  study 
the  diagnostic  fields  for  varying  initial  times.  The  structures 
observed  in  these  fields  are  related  to  different  mixing  pro¬ 
cesses  occurring  during  the  three  main  phases  of  the  model: 
crystalization  in  which  polygonal  eyewall  features  form  and 
develop  filamentation,  vortex  interaction  and  merger  which 
destroy  the  symmertry,  and  final  collapse  into  a  monopole. 
Section  6  is  devoted  to  the  mixing  rates  mentioned  before, 
which  are  displayed  as  functions  of  initial  time  and  either  ra¬ 
dius  or  value  of  S  along  a  shear-line.  Concluding  remarks 
and  an  outlook  on  future  work  are  given  in  Sect.  7. 


2  Model  overview 


The  model  used  in  this  paper  is  the  2-D  nondiver¬ 
gent  barotropic  model  for  hurricane-like  vortices  studied 
by  Kossin  and  Schubert  (2001);  Kossin  and  Eastin 
(2001);  Schubert  et  al.  (1999).  The  velocity  field, 
u (x,t)=(u(x,t),v(x,t))*  with  x=(x,y)*  eM2  (asteriks  de¬ 
note  transposed  vectors  or  matrices),  is  given  as  the  solution 
on  the  /-plane  of  the  incompressible  Navier  Stokes  equation, 

—  +  (u-V)u-/Bu+!vf>  =  vV2u,  (1) 

3 1  p 

V  u  =  0,  (2) 

where 


p  is  the  pressure,  p  the  constant  density,  /  the  constant  Cori¬ 
olis  parameter,  and  v  the  constant  viscosity,  chosen  to  be 
100  m2  s_1.  In  the  choice  of  v  we  follow  Schubert  et  al. 
(1999),  while  Kossin  and  Schubert  (2001)  used  5m2s_1. 
The  choice  of  viscosity  may  have  an  effect  on  long  time 
mixing  processes,  which  could  be  studied  by  the  methods 
of  Hendricks  and  Schubert  (2009).  Expressing  the  velocity 
in  terms  of  a  streamfunction  as  u=— BV^r  and  elim¬ 

inating  the  pressure  from  Eq.  (1),  leads  to  the  equation 


3£  dx/r  3£  dx/r  3£ 
3 1  dx  3 y  3 y  dx 


(3) 


where  ^=V2^  is  the  relative  vorticity.  Following  Kossin  and 
Schubert  (2001),  we  impose  periodic  boundary  conditions 
on  xjs  with  a  fundamental  domain  of  600  km  x  600  km,  and 
choose  as  initial  condition  an  almost  circular  symmetric  ring 
of  vorticity,  ^o(r,0),  to  model  a  2-D  hurricane  after  an  initial 
eyewall  has  formed.  The  defining  equation  of  £o (r,6)  is  the 
equation  used  in  Kossin  and  Schubert  (2001). 

Equation  (3)  was  solved  numerically  using  a  Fourier 
pseudospectral  method  with  512x512  collocation  points. 
Dealiasing  results  in  170x170  Fourier  modes.  The  ODE- 
system  for  the  Fourier  modes  was  solved  via  Matlab’s 
ode45  routine,  which  implements  a  fourth  order  Runge-Kutta 
method  with  adaptive  time  steps. 

In  our  numerical  calculation  of  x/r  and  f ,  we  reproduced 
the  behavior  observed  in  Schubert  et  al.  (1999).  The  annular 
ring  of  high  vorticity  fluid  develops  a  wavenumber  4  asym¬ 
metry,  which  is  present  in  the  vorticity  fields  as  early  as  2  h, 
and  develops  into  a  polygonal  eyewall,  with  4  mesovortices 
after  6  h.  After  8  h,  the  mesovortices  begin  to  break  down  and 
merge.  The  breakdown  of  the  mesovortices  is  nearly  com¬ 
plete  after  12  h,  and  mixing  of  high  and  low  vorticity  occurs 
along  long  filament  structures.  The  relative  vorticity  fields 
during  these  times  can  be  seen  in  Figs.  6  to  11a.  After  24  h, 
diffusive  mixing  along  the  filaments  leads  to  a  more  mixed 
state.  Few  pools  of  high  or  low  vorticity  fluid  remain,  with  a 
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pool  of  low- vorticity  fluid  from  the  eye  migrating  through  the 
eyewall,  and  high  vorticity  fluid  redistributing  in  the  eye.  Af¬ 
ter  48  h,  a  high  vorticity  eye  and  a  low  vorticity  environment 
remain  in  a  monopole  endstate.  The  eyewall  is  no  longer 
present  as  there  is  no  longer  a  strong  angular  velocity  gradi¬ 
ent. 

In  this  paper  we  study  Lagrangian  mixing  in  the  model, 
which  is  based  on  following  trajectories  for  varying  initial 
times.  The  trajectories  were  calculated  with  the  same  spatial 
and  temporal  resolution  as  the  model  output,  using  a  fourth 
order  Runge-Kutta  method  with  a  fixed  time  step  of  7.5  s. 
Because  of  time  and  memory  limitations  associated  with  the 
large  number  of  trajectories  needed  for  quantifying  mixing 
over  a  sequence  of  initial  times,  the  trajectories  used  for  com¬ 
puting  time-dependent  mixing  rates  were  calculated  with  a 
time  step  of  60  s.  Comparison  of  the  results  for  the  two  time 
steps  for  a  small  random  set  of  initial  conditions  showed  that 
the  use  of  the  coarser  time  resolution  in  the  mixing  calcula¬ 
tions  is  justified. 

3  Diagnostic  fields  for  mixing  and  shear 

In  this  section  we  introduce  the  scalar  fields  utilized  to  di¬ 
agnose  the  particle  flow  resulting  from  the  numerically  cal¬ 
culated  velocity  field.  A  main  characteristic  feature  of  the 
model  is  an  almost  circular  motion,  the  trajectories  encir¬ 
cle  the  origin  in  the  counterclockwise  direction.  The  model 
shows  a  strong  variation  of  the  particle  speed  |u|  in  the  radial 
direction.  This  variation  leads  to  high  shearing  that  dom¬ 
inates  the  particle  separation,  but  is  not  the  result  of  hyper¬ 
bolic  trajectory  separation.  Superimposed  on  this  shear  effect 
is  hyperbolic  mixing  due  to  trajectory  separation  in  directions 
orthogonal  to  the  velocity. 

In  order  to  diagnose  hyperbolicity,  we  exploit  the  La¬ 
grangian  field  introduced  in  Haller  and  Iacono  (2003),  in 
which  hyperbolic  trajectory  splitting  is  separated  from  par¬ 
ticle  separation  due  to  shearing.  The  more  common  FTLE- 
field  is  also  analyzed,  however,  this  field  is  dominated  by  the 
shear  and  hence  not  suitable  as  an  indicator  for  hyperbolic 
mixing.  In  order  to  quantify  hyperbolic  mixing,  we  define 
closed  shear  lines  as  contour  lines  of  a  suitably  defined  shear 
field,  and  measure  and  predict  transport  across  these  lines 
(Sect.  6).  Further  indicators  used  in  our  study  are  two  Eu- 
lerian  fields:  The  Hessean  determinant  of  the  streamfunction 
(2-field),  and  the  radial  gradient  of  the  instantaneous  angular 
velocity. 

3.1  Eulerian  fields 

3.1.1  Q-field 

Eulerian  trajectory  separation  occurs  when  the  linearized  ve¬ 
locity  shows  local  expansion  of  area.  The  local  variation  of 
area  can  be  inferred  from  the  Jacobian  of  the  velocity  field, 


that  is,  from  the  Hessian  determinant  of  the  streamfunction, 
which  is  referred  to  as  the  Q -field, 

Q(X,t)  =  ^y(X,t)  -  fXX(X,t)l/fyy(X,t).  (4) 

According  to  the  Okubo- Weiss  criterion  (Schubert  et  al., 
1999),  regions  with  2>0  show  local  trajectory  repulsion, 
whereas  regions  with  g<0  show  local  attraction.  The  Q- 
field  allows  diagnosis  of  instantaneous  separation,  which  typ¬ 
ically  differs  from  Lagrangian  measures  of  separation. 

3.1.2  Angular  velocity 

The  strong  rotation  and  near  symmetry  of  the  flow  suggests 
that  polar  coordinates  (r,0)  provide  a  useful  coordinate  sys¬ 
tem  for  displaying  fields  calculated  from  the  velocity  field.  In 
particular,  the  quasi-circular  behavior  of  trajectories  suggests 
that  the  angular  velocity,  <z>=r_1u-(—  sin0,cos0)*,  is  an  ap¬ 
proximate  measure  of  the  particle  speed,  and  the  derivative 
dco/dr  is  an  approximate  measure  of  shearing. 

For  any  scalar  field  (p(x,  t ),  a  measure  for  the  radial  varia¬ 
tion  is  provided  by  the  angular  average,  indicated  by  an  over¬ 
bar, 

i  r27T 

(p(rj)=—  /  <p(r,6,t)d6 . 

2tt  Jo 

Contours  of  are  shown  in  Figs.  6  to  1  ld-f  showing  the  re¬ 
lationship  of  maxima  (maximum  normal  propogating  shear), 
and  minima  (maximum  counter  propogating  shear)  of  to 
features  of  other  scalar  fields. 

3.2  Lagrangian  fields 

Let  (xo)  be  the  flow  map  associated  with  the  equation 

x  =  u(x,0  (5) 

for  particle  trajectories  x(t),  that  is,  the  solution  of  Eq.  (5) 
with  initial  condition  x(to)=xo.  Small  perturbations  in  the 
initial  condition,  y0=xo+§o>  lead  to  a  perturbed  trajectory 
y (t)=x(t)+i;(t).  For  sufficiently  small  |§ol ,  the  perturbation 
§(0  can  be  approximated  through  the  Jacobian  of  the  flow 
map  as 

m  =  V<(xo)!o,  (6) 

which  satisfies  the  variational  equation 
$  =  Vu(x(f),0$.  (7) 

3.2.1  Finite  time  Lyapunov  exponents 

Finite  time  Lyapunov  exponents  have  become  a  standard  in¬ 
dicator  for  Lagrangian  trajectory  separation  (Haller  and  Poje, 
1997;  Haller,  2000;  Haller  and  Yuan,  2000;  Haller,  2002). 
Consider  a  time  range  [to,to+T]  with  fixed  integration  time 
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T.  The  growth  of  |£(f)|  during  this  time  range  is  governed 
by  the  Cauchy-Riemann  deformation  tensor, 

A(xo,to)=(dxo0!o+r(xo))*(dxoK+T(x°))’ 


fluid  velocity  is  non- autonomous,  its  time  variation  is  slow, 
so  that  u (x(t),t)  is  still  close  to  a  solution  of  Eq.  (7)  for  finite 
times.  The  transformed  system  for  rj  can  be  written  in  the 
form  (Haller  and  Iacono,  2003), 


and  becomes  maximal  when  §o  is  aligned  with  the  eigenvec¬ 
tor  corresponding  to  the  larger  eigenvalue,  Amax(A),  of  A. 
The  quantity 

1 

a  (x0 ,  *o) = In  kmax  (  A  (x0 ,  *o)) 

is  the  FTLE -field  (considered  as  a  function  of  xo)  at  the  ini¬ 
tial  time  to.  One  distinguishes  forward  TTLE-fields  (T >0) 
and  backward  FTZJs-fields  (T < 0).  Maximal  ridges  of  the 
forward  and  backward  TTLE-fields  play  the  role  of  unsta¬ 
ble  (repelling)  and  stable  (attracting)  manifolds  over  the  time 
range  considered  (Haller,  2000;  Shadden  et  al.,  2005).  In  the 
limit  T^oo  or  T^—o o,  a  approaches  the  true  Lyapunov 
exponent,  and  the  maximal  ridges  approach  the  stable  and 
unstable  manifolds,  of  the  full  trajectory  seeded  at  (xo,fo)- 
Data  limitations  in  the  case  of  numerically  computed  veloc¬ 
ity  fields  allow  only  finite  time  integrations,  and  restrict  the 
analysis  to  FTLE’s. 

For  our  velocity  field,  the  particle  separation  is  dominated 
by  the  shearing  in  the  radial  direction.  As  a  result,  the  FTLE- 
values  measure  growth  of  perturbations  in  approximately  an¬ 
gular  directions,  and  high  FTLE- values  (ridges)  occur  near 
extrema  of  dco/dr ,  whereas  low  FTLE- values  occur  near 
zero  contours  of  dco/dr.  Generally,  FTLE- fields  are  not  suit¬ 
able  as  indicators  of  hyperbolic  mixing  in  the  presence  of 
high  shear. 

3.2.2  Integrated  Q -field 


?  =  [A(0  +  K0B>, 

where 


(10) 


A(0= 


V  0 


t)  a(x(t),t )" 

r(x(0,0  7’ 


1 


*(0  =  71 2U  *ulx=x(0> 


r(x,0=T — ^(ui)*(Vu)u1 


=  i(u2Vy-UV(Uy  +  VX)  +  V2UX), 

\u\z  \  / 

and  the  non-diagonal  entry  a  is  composed  of  two  parts, 
a(x,t)=s(x,t)+d(x,t ), 
where 

s(x,0=2-u*(Vu)u-l 


=  ■— 2  (u2Uy+UV(Vy—Ux)  —  V2Vx), 


d(x,t)=— ,(u7(Vu)u 


=  — j(u2VX+UV{Vy-UX)-V2Uy). 


In  addition  to  the  instantaneous  Q -field,  Eq.  (4),  we  consider 
the  integrated  Q -field,  formed  by  integrating  Q  along  trajec¬ 
tories 

^  fto+T 

G(x,0=  /  Q(xo(t),t)dt.  (8) 

Jt0 

3.2.3  Lagrangian  fields  for  hyperbolic  mixing  and  shear 

Following  Haller  and  Iacono  (2003),  in  order  to  separate 
mixing  and  shear  in  the  variational  system  (Eq.  7),  a  mov¬ 
ing  frame  of  reference  is  introduced  by  setting 

I  =  M(x(t),t)ri,  (9) 

where  the  component  vectors  of  the  matrix  M, 

M(x,0=  1  (u(x,r),uJ-(x,r)),  u_L=(— 

|u(x,0I 

are  the  normalized  fluid  velocity  u/|u|,  and  the  unit  vector 
orthogonal  to  u.  This  transformation  is  motivated  by  the 
fact  that  for  autonomous  velocity  fields,  u=u(x),  u (x(t))  is 
a  solution  of  Eq.  (7).  Although  our  numerically  computed 


The  terms  in  the  transformed  linearized  system  (Eq.  10)  mo¬ 
tivate  the  definition  of  Lagrangian  fields  as  diagnostics  for 
hyperbolic  mixing  and  shear.  Since  our  velocity  field  is 
slowly  varying  in  time,  the  terms  associated  with  b(t)  in 
Eq.  (10)  are  neglected  in  these  definitions. 

3.3  R-field 

As  a  consequence  of  incompressibility,  the  matrix  A (t)  has 
the  eigenvalues  ±r(x(f),  t).  Fixing  an  integration  time  T,  the 
integrated  field  R , 

R(xo,to)=  /  r(<pt°+T (x0),t)  dr,  (11) 

describes  the  growth  of  a  perturbation  in  the  direction  orthog¬ 
onal  to  the  Lagrangian  velocity,  and  the  ratio  R/\T\  plays 
the  role  of  a  finite-time  Lyapunov  exponent  in  this  direction. 
Thus  R  is  a  measure  of  attraction  (R<0)  or  repulsion  (R>0) 
of  nearby  trajectories  towards  x(t)  over  the  integration  inter¬ 
val  [to,to+T].  Due  to  incompressibility,  expansion  orthog¬ 
onal  to  u  is  combined  with  contraction  in  the  direction  of  u 
and  vice  versa. 
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We  note  that  in  the  limit  T^oo,  R  reduces  to  the  mixing 
efficiency  proposed  by  Ottino  (1989),  when  this  efficiency  is 
evaluated  in  the  direction  orthogonal  to  u.  In  our  study,  R 
will  be  used  as  the  main  diagnostic  field  for  hyperbolicity, 
and  in  addition  as  a  means  to  predict  mixing  rates  across  the 
shear-lines  defined  below. 

3.4  S-field 

The  term  s(x,  t)  can  be  written  in  the  form 

s(x,t)  -  (V|u(x,r)|)-u-L(x,0,  (12) 

and  hence  characterizes  the  rate  of  change  of  the  particle 
speed  in  the  direction  orthogonal  to  the  velocity.  Thus  s  is 
a  local,  Eulerian  measure  of  shear  in  the  fluid  flow.  We  de¬ 
fine  the  S-field  by  integrating  s  along  trajectories, 

Ct0  +  T  t  JL 

S(x0,to)=  s((pt°+r(x0),T)  dr,  (13) 

Jt0 

and  use  S  as  a  Lagrangian  diagnostic  field  for  shearing. 

We  note  that  an  alternative  Lagrangian  measure  of  shear 
has  been  defined  in  Haller  and  Iacono  (2003)  using  the  non¬ 
diagonal  entry  of  the  fundamental  matrix  of  ij  =  A (t)rj. 
However,  this  field  involves  a  double  time-integral  and  is 
computationally  more  expensive.  The  S-field  has  a  straight¬ 
forward  interpretation  as  shear-diagnostic  due  to  Eq.  (12), 
and  requires  less  computational  effort. 

As  in  the  case  of  the  TTLF-field,  we  distinguish  forward 
(r>0)  and  backward  (T<0)  fields  for  both  R  and  S. 

3.5  Shear-lines 

For  a  given  integration  time  T,  we  define  the  shear- line  of 
strength  C  at  initial  time  to  as  the  level  contour  of  S,  i.e., 

5,c={xo|5,(xo,ro)=C'} 

.  High  values  of  |C|  correspond  to  lines  with  high  shear. 
For  our  model,  the  shear  lines  are  all  closed  curves  around 
the  origin  (distorted  circles).  Positive  and  negative  values  of 
S  indicate  that  the  Lagrangian  speed  increases  when  mov¬ 
ing  radially  outwards  and  inwards,  respectively.  We  refer  to 
the  first  case  as  “normal  propagating  shear”  and  to  the  sec¬ 
ond  case  as  “counter-propagating  shear”.  Hyperbolic  mixing 
measured  by  R  is  associated  with  transport  across  the  shear¬ 
lines.  This  will  be  used  in  Sect.  6  to  define  mixing  rates. 


4  Lagrangian  fields  and  coherent  structures 

In  Figs.  1,  2,  and  3  we  show  forward  and  backward  FTLE-, 
S- ,  and  R-fields,  respectively,  at  the  initial  time  ^o=6h,  af¬ 
ter  the  polygonal  eyewall  has  formed,  and  for  integration 
times  r=zbl5min,  r=ib30min,  and  T= ibl20min.  For 
r=zbl5min,  the  FTTF’-field  (Fig.  la,b)  reveals  coherent 


structures  near  the  polygonal  eyewall,  since  the  effect  of  the 
shear  is  not  so  pronounced  over  this  short  time  range.  For 
increasing  |  T  |  the  shear  becomes  dominant,  and  the  FTLE- 
level  contours  evolve  into  distorted  circles  (Fig.  lc-f).  Com¬ 
parison  of  the  azimuthally  averaged  fields  (Sect.  5)  shows 
that  high  FTLE- values  occur  near  extreme  values  of  dco/dr. 
While  high  FTLE- values  correspond  to  high  trajectory  sepa¬ 
ration,  they  do  not  give  clear  LCS’s  (ridges)  for  longer  inte¬ 
gration  times.  LCS’s  can  be  seen  only  at  very  short  integra¬ 
tion  times  (Fig.  la  and  b),  and  in  regions  that  are  predicted 
by  the  Q -field.  As  integration  time  is  increased  (Fig.  lc- 
f),  the  LCS’s  do  lengthen  as  expected,  but  they  also  become 
broader.  In  particular,  the  LCS’s  for  T— ±  15  min  that  are  lo¬ 
cated  near  the  corners  of  the  eyewall,  converge  into  a  single 
broad  ring  that  represents  an  annulus  of  high  shear. 

The  S-field  (Fig.  2)  is  a  shear-indicator,  and  its  level  con¬ 
tours  (the  shear-lines)  are  distorted  circles  for  all  integration 
times.  The  TTLF-field  shows  similar  structures  as  the  S-field 
for  longer  integration  times,  confirming  that  trajectory  sepa¬ 
ration  is  mainly  due  to  shear. 

The  R -field  (Fig.  3)  shows  structures  of  high  and  low  R- 
values  that  persist  over  a  series  of  integration  times,  making 
them  coherent.  These  structures  lengthen  and  become  more 
resolved  (narrower)  when  the  integration  time  increases.  Ini¬ 
tial  points  on  ridges  and  valleys  have  R> 0  and  R< 0,  indicat¬ 
ing  strong  separation  and  contraction  in  the  (approximately 
radial)  direction  orthogonal  to  the  Lagrangian  velocity,  re¬ 
spectively.  The  structures  exist  in  both  the  forward  and  back¬ 
ward  time  fields,  and  some  of  the  forward  and  backward  time 
structures  have  intersection  points.  Since  the  R -field  is  radi¬ 
ally  continuous,  high  values  of  the  R-field  lead  to  trajectories 
that  show  high  net  movement  orthogonal  to  the  Lagrangian 
velocity,  and  hence  are  more  likely  to  cross  shear-lines.  Since 
the  structures  span  across  the  shear-lines,  they  are  not  La¬ 
grangian,  as  trajectories  with  high  angular  velocity  pass  tra¬ 
jectories  with  lower  angular  velocity. 

A  prominent  feature  of  the  forward  R -field  at  to=6h  (and 
later)  are  the  filaments  observable  in  Fig.  3a,  c  and  e,  which 
are  a  consequence  of  the  polygonal  eyewall.  At  ^o=2h 
(Fig.  4a)  no  filaments  are  observed.  The  filamentation  con¬ 
cerns  the  ridges  and  valleys,  as  well  as  the  edges  between 
them. 

4.1  An  advective  mixing  mechanism 

Initial  conditions  xo  that  satisfy  /?(xo,/o)=0  are  invariant  in 
the  sense  that  there  is  no  net  movement  of  neighboring  tra¬ 
jectories  relative  to  the  Lagrangian  velocity.  As  can  be  seen 
in  Figs.  3  and  4a,  ridges  and  valleys  of  R  come  in  nearby 
pairs,  and  the  ridge  and  valley  of  a  pair  are  separated  by  a 
segment  of  a  zero  contour  which  forms  (approximately)  an 
edge  of  the  R -field.  The  edge  is  neutrally  stable,  that  is,  it 
attracts  from  one  side  (from  the  ridge)  and  repels  from  the 
other  side  (towards  the  valley).  The  situation  is  illustrated  in 
Fig.  4a  ,c  and  d  for  to=2h.  At  this  initial  time  the  motion  is 
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Fig.  1.  Forward  and  backward  FTLE-fields  at  the  initial  time  t$  =  6  hours  integrated  (a)  15  min,  (b)  —15  min,  (c)  30  min,  (d)  —30  min,  (e) 
120  min,  and  (f)  —120  min. 
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Fig.  2.  Forward  and  backward  5-fields  at  the  initial  time  t$=6\\  integrated  (a)  15 min,  (b)  —15 min,  (c)  30 min,  (d)  —30 min,  (e)  120 min, 
and  (f)  —120  min. 
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Fig.  3.  Forward  and  backward  R-fields  at  the  initial  time  to=6h  integrated  (a)  15  min,  (b)  —15  min,  (c)  30  min,  (d)  —30  min,  (e)  120  min, 
and  (f)  —120  min. 
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Fig.  4.  (a)  R- field  with  ridges  and  valleys  and  (b)  Lagrangian  speed  at  t$=2h.  (c)  Zoom  of  box  A  from  /?-field  showing  repelling  edge,  and 
(d)  zoom  of  box  B  showing  attracting  edge.  Integration  time  is  T =60  min.  Black  lines  show  the  azimuthal  velocity  at  the  initial  time. 


almost  circular  (Fig.  4b),  but  motion  across  shear-lines  can 
be  observed  already. 

Let  TZ  be  a  structure  (ridge  or  valley)  of  R  at  initial  time 
to.  The  structure  is  coherent  in  the  sense  that  it  evolves  con¬ 
tinuously,  for  varying  initial  time,  into  a  structure  TZ'  seen  as 
a  ridge  or  valley  of  R  at  initial  time  to+x.  It  is,  however, 
not  Lagrangian  because  it  is  not  advected  with  the  flow,  that 
is,  the  image  of  7 Z  under  the  flow  map,  7£r=0^+r  (TZ),  has 
advanced  farther  from  7 Z  than  TZ',  and  is  not  a  structure  of 
R  at  to+T  (Fig.  5a).  Generally  we  observe  that  the  coherent 
structures  move  at  a  slower  rotational  speed  than  that  of  the 


mean  flow,  which  can  be  attributed  to  the  effect  of  Rossby 
waves  (Montgomery  and  Lu,  1997). 

A  coherent  structure  TZ  has  a  leading  and  a  lagging  end  rel¬ 
ative  to  counterclockwise  rotation.  Concerning  the  evolution 
of  TZ  under  the  flow  map,  two  cases  can  occur  for  structures 
computed  in  a  forward  time  integration: 

(a)  If  the  leading  end  is  at  higher  angular  velocity  than  the 
lagging  end,  then  the  image  TZj  of  TZ  under  0^+r  is 
lengthened  over  the  integration  and  tends  to  align  with 
a  contour  of  the  S-field  (Fig.  5b). 
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(b)  If  the  leading  end  is  at  lower  angular  velocity  than  the 
lagging  end,  then  the  flow  map  rotates  7 2,  and  for  suf¬ 
ficiently  large  T  the  image  7 2t  tends  to  align  with  a 
contour  of  the  S-field  in  the  opposite  direction  (Fig.  5c). 

For  a  nearby  pair  of  a  ridge  and  a  valley,  the  relative  posi¬ 
tion  of  their  flow  map  images  is  preserved  in  case  a,  whereas 
in  case  b  they  switch  position.  This  rotation  and  position 
switching  are  a  mechanism  for  the  advective  mixing  during 
the  polygonal  eyewall  stage.  The  square  eyewall  gives  four 
valleys  and  ridges  aligned  in  a  way  that  four  sets  of  trajec¬ 
tories  pass  from  outside  to  inside  and  four  from  inside  to 
outside  of  an  ^-contour.  The  combination  of  ridge-edge- 
valley  sets  of  R ,  aligned  with  the  leading  ends  at  lower  angu¬ 
lar  velocity,  can  be  seen  as  an  indicator  of  fluid  regions  that 
will  roll  into  mesovortices  over  the  forward  integration  time. 
Thus  the  strength  and  size  of  the  surrounding  ridges  and  val¬ 
leys  are  an  indicator  of  the  potential  flux  in  and  out  of  the 
mesovortex. 

5  Field  diagnostics  for  varying  initial  time 

In  Figs.  6  to  11  we  show  in  panel  a  the  relative  vorticity,  in 
panel  b  the  Q -field,  and  in  panel  c  the  R -field  with  over- 
layed  vorticity  contours  chosen  to  illustrate  the  relation  be¬ 
tween  R -field  structures,  and  vorticity  structures.  We  show 
the  S-field,  integrated  Q -field,  and  FTLE  field  in  panel  d, 
e,  and  f  respectively,  together  with  contours  (from  inside 
to  outside)  of  the  maximum  normal  propogating  shear,  the 
maximum  tangential  velocity,  and  the  maximum  counterpro- 
pogating  shear.  The  initial  times  in  these  figures  are  to=2  h, 
4  h,  6  h,  8  h,  10  h,  and  12  h,  and  the  integration  time  is  T=1  h. 

Since  the  shear-lines  are  distorted  circles,  we  can  interpret 
the  average  a(ro,to)  (ro=|xo|)  as  radial  mixing  rate.  In  all 
Figs.  6  to  1  Id— f  we  observe  that  extreme  Q ,  S ,  and  FTLE- 
values  occur  at  extrema  of  dco/dr,  demonstrating  that  the  Q , 
S ,  and  FTLE -fields  are  dominated  by  the  shear.  A  similar 
interpretation  as  radial  mixing  rates  can  be  attributed  to  the 
averages  S(ro,  to)  and  R(ro,  to).  Plots  of  S  and  a  reveal  these 
two  averages  are  very  similar  in  structure,  as  both  measure 
shear.  The  quantity  R  can  be  interpreted  as  a  measure  of 
hyperbolic  mixing,  which  is  important  for  transport  through 
the  eyewall. 

5.1  2-4  h:  initial  state 

At  the  initial  time  of  2h  (Fig.  6),  the  model  is  still  close 
to  the  initial  state  and  shows  a  broad  ring  of  high  vorticity 
fluid.  While  the  vorticity,  Q,  S,  and  FTLE- fields  are  almost 
circular-symmetric,  the  R-field  shows  distinct  lines  of  high 
radial  mixing,  demostrating  that  coherent  structures  can  per¬ 
sist  through  dominant  shear.  The  wavenumber  four  asym¬ 
metry  begins  to  show  in  the  R -field,  particularly  in  the  for¬ 
ward  time  integration,  atthough  the  initial  vorticity  profile  is 
nearly  preserved,  with  any  asymmetries  barely  noticible.  The 


i  I 


S  contour 


(a) 


(c) 


Fig.  5.  (a)  Sketch  of  a  structure  12  of  R  at  initial  time  to ,  the  associ¬ 
ated  structure  12r  of  R  at  initial  time  to+r,  and  the  flow  map  image 
72r=0^+T(72),  illustrating  the  non-Lagrangian  nature  of  the  co¬ 
herent  structures,  (b)  Structure  12  crossing  an  S  contour  with  speed 
above  the  contour  higher  than  the  speed  below,  and  flow  map  image 
12j  after  an  integration  time.  Leading  and  lagging  ends  of  12  and 
their  images  on  12  j  are  marked  by  a  triangle  and  a  square,  respec¬ 
tively.  (c)  Same  as  (b)  with  opposite  orientation  of  12  relative  to  the 
contour,  leading  to  a  rotation  of  12j . 


www.  atmos-chem-phy  s  .net/ 1 0/47 5/20 1 0/ 


Atmos.  Chem.  Phys.,  10,  475-497,  2010 


486 


B.  Rutherford  et  al.:  Advective  mixing  in  a  nondivergent  barotropic  hurricane  model 


Fig.  6.  (a)  Relative  vorticity  field,  (b)  Q -field,  (c)  R -field  with  vorticity  contours  overlayed  to  relate  structures,  (d)  S-field,  (e)  integrated 
2-field  Q,  and  (f)  TTLF-field  for  initial  time  of  tQ=2h  with  integration  time  r=lh.  The  white  contours  in  (d)-(f)  mark  the  radius  of 
maximum  (inner,  solid),  the  radius  of  maximum  tangential  wind  (middle,  dashed),  and  the  radius  of  maximum  counterpropogating 
(outer,  dash-dot) 
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rotation  of  the  R -ridges  and  valleys  in  this  stage  allows  the 
crystallization  that  is  neccessary  for  mesovortex  formation. 

At  4  h  (Fig.  7),  much  of  the  symmetry  of  the  initial  state 
still  remains.  A  forward  time  integration  begins  to  display 
asymmetries  in  the  Q ,  FTLE-  and  S-fields,  whereas  the  R- 
field  retains  the  structures  that  were  present  at  2h.  At  4h 
the  backward  time  integration  of  R  (not  shown)  also  shows 
regions  of  high  orthogonal  separation. 

5.2  6-8  h:  polygonal  eye  wall 

At  6  h  (Fig.  8),  the  vorticity  field  shows  a  polygonal  eyewall 
structure,  where  the  flow  resembles  a  nonlinear  critical  layer 
for  dry  barotropic  instability.  A  square  inner  eyewall  struc¬ 
ture  forms,  with  pools  of  low  vorticity  fluid  organized  into 
the  corners  of  the  eye.  The  fluid  is  largely  unmixed,  with 
low  vorticity  fluid  organized  in  the  mesovortices,  where  it  is 
largely  protected  from  the  outer  flow.  Since  vorticity  is  ma¬ 
terially  conserved,  low  vorticity  fluid  from  the  eye  and  en¬ 
vironment  must  replace  the  fluid  that  left  the  eyewall.  The 
fluid  that  is  mixed  across  the  boundaries  is  consistent  with 
the  ridges  and  valleys  of  the  R -field.  The  filamentation  that 
develops  from  the  stretching  of  high  vorticity  fluid  that  exits 
from  the  vorticity  ring  to  the  environment,  can  be  seen  in  the 
form  of  spiral  bands  in  the  vorticity  field. 

The  Q -field  shows  that,  instantaneously,  high  trajectory 
separation  occurs  along  the  square  boundary  of  the  inner  eye- 
wall.  The  pools  of  low  vorticity  show  instantaneous  con¬ 
traction.  The  outer  vorticity  ring  shows  high  trajectory  sep¬ 
aration.  Even  when  the  FTLF-fields  are  calculated  for  the 
small  integration  time  of  T =3  min,  there  is  a  noticible  dif¬ 
ference  between  the  separation  points  of  the  Q -field  and  the 
FTLF-field.  The  square  eyewall  formation  corresponds  to 
four  structures  of  high  FTLE- values  in  both  the  forward  and 
backward  time  fields  for  short  integration  times.  As  inte¬ 
gration  time  is  increased,  the  structures  lengthen  and  are  no 
longer  distinguishable. 

The  R -field  shows  a  series  of  ridges  and  valleys  that  orig¬ 
inated  as  coherent  structures  from  the  earlier  times,  but  are 
not  as  refined  as  previous  structures.  There  are  also  structures 
emanating  outward  from  the  ring  of  high  vorticity  that  may 
play  the  role  of  protecting  the  ring  from  interaction  with  the 
outer  flow  (Dunkerton  et  al.,  2009).  If  our  model  was  a  true 
representation  of  a  wave  critical  layer,  the  structures  would 
correspond  to  dividing  streamlines,  but  the  high  shear  pro¬ 
hibits  this.  The  maximal  R -regions  are  located  at  the  same 
places  that  show  high  FTLE- values,  but  with  much  greater 
resolution  than  the  FTLF-fields,  which  are  blurred  by  the 
shear. 

At  8h  (Fig.  9),  the  polygonal  eyewall  structure  that  is 
present  at  6  h  is  still  clearly  visible.  The  R-field  has  simi¬ 
lar  properties  as  the  R-field  at  6h  (Fig.  8),  with  the  coherent 
structures  begining  to  merge,  showing  intense  mixing. 


5.3  Later  state:  mixing  into  a  monopole 

The  period  after  to=8  h  until  ^o=12h  exhibits  intense  mixing 
that  leads  to  a  collapse  into  a  monopole  end  state. 

For  to=10h  (Fig.  10),  the  g-field  is  less  square,  and  the 
R -field  shows  high  mixing  in  two  distinct  regions,  one  with 
expansion  and  one  with  contraction.  The  region  of  expan¬ 
sion  is  inside  the  ring  of  high  vorticity,  while  the  region  of 
contraction  has  become  organized  in  the  dominant  mesovor¬ 
tex,  which  is  the  “winner”  and  survives  to  become  the  pri¬ 
mary  vortex  during  the  collapse  into  the  monopole  end  state. 
Note  that  the  merging  of  LCS’s  into  a  single  LCS  is  a  bifur¬ 
cation,  and  cannot  happen  if  they  represent  true  stable  and 
unstable  manifolds,  however  the  coherent  structures  here  are 
not  entirely  Lagrangian,  yet  their  interactions  and  bifurca¬ 
tions  play  an  important  role  in  the  systematic  mixing  during 
mesovortex  interaction.  As  more  advection  (stirring)  in  and 
out  of  the  eyewall  occurs,  there  is  filamentation  of  the  initial 
vorticity-contours  with  diffusive  mixing  occurring  along  the 
lengthening  contour  boundaries,  leading  to  an  “averaging”  of 
vorticity  values  through  diffusion. 

For  to=12h  (Fig.  11),  the  inner  ring  of  vorticity  has  bro¬ 
ken  down.  The  FTLE-  and  ^-fields  show  the  outer  rings 
converged  as  a  thick  ring,  and  the  model  is  entering  the 
monopole  state.  The  R -field  coherent  structures  now  show 
the  dominant  mesovortex  migrating  to  the  center  and  the 
other  mesovortices  are  disappearing  due  to  their  annihilation 
by  the  dominant  mesovortex.  Regions  of  high  R -values  are 
pushed  outward,  indicationg  mixing  with  the  outer  flow.  At 
this  stage,  there  is  a  single  remaining  protecting  R -ridge,  on 
the  outside  of  the  remaining  high  vorticity  ring,  which  has 
served  the  role  of  protecting  the  mesovortex  that  eventually 
becomes  the  “winner”. 

Beyond  12  h,  the  initial  regions  of  vorticity  are  not  recog¬ 
nizable,  and  high  (although  not  as  high  as  the  initial  state) 
vorticity  fluid  begins  to  organize  into  the  eye.  The  low  vor¬ 
ticity  fluid  from  the  eye  becomes  well  mixed,  and  the  eyewall 
and  environment  become  filled  with  relatively  low  vorticity 
fluid.  The  angular  velocity  gradient  decreases,  and  the  S- 
field  shows  no  eyewall. 

Although  the  fluid  that  is  mixed  beyond  12  h  is  not  distin¬ 
guishable  based  on  its  initial  vorticity,  the  R-field  still  gives 
regions  of  advective  mixing,  showing  that  the  moving  frame 
of  initial  conditions  still  shows  regions  of  fluid  that  are  trans¬ 
ported. 


6  Mixing  rates 

The  radial  mixing  rates  a(ro,^o)  and  S(roTo)  quantify  mix¬ 
ing  due  to  shear,  whereas  R(ro,to)  quantifies  hyperbolic  mix¬ 
ing.  In  these  mixing  rates,  the  lines  along  and  across  which 
mixing  is  quantified  are  circles.  Hyperbolic  mixing  rates  that 
are  more  closely  related  to  the  shearing  structures  are  mixing 
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Fig.  7.  (a)  Relative  vorticity  field,  (b)  Q -field,  (c)  R -field  with  vorticity  contours  overlayed  to  relate  structures,  (d)  S-field,  (e)  integrated 
2-field  Q ,  and  (f)  TTLE-field  for  initial  time  of  tQ=4h  with  integration  time  r=lh.  The  white  contours  in  (d)-(f)  mark  the  radius  of 
maximum  (inner,  solid),  the  radius  of  maximum  tangential  wind  (middle,  dashed),  and  the  radius  of  maximum  counterpropogating 
(outer,  dash-dot). 
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Fig.  8.  (a)  Relative  vorticity  field,  (b)  Q -field,  (c)  R -field  with  vorticity  contours  overlayed  to  relate  structures,  (d)  S-field,  (e)  integrated 
2-field  Q,  and  (f)  TTLF-field  for  initial  time  of  tQ=6  h  with  integration  time  r=lh.  The  white  contours  in  (d)-(f)  mark  the  radius  of 
maximum  (inner,  solid),  the  radius  of  maximum  tangential  wind  (middle,  dashed),  and  the  radius  of  maximum  counterpropogating 
(outer,  dash-dot). 
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Fig.  9.  (a)  Relative  vorticity  field,  (b)  Q -field,  (c)  R -field  with  vorticity  contours  overlayed  to  relate  structures,  (d)  S-field,  (e)  integrated 
g-field  Q ,  and  (f)  FTLE- field  for  initial  time  of  fy-Sh  with  integration  time  r=lh.  The  white  contours  in  (d)-(f)  mark  the  radius  of 
maximum  (inner,  solid),  the  radius  of  maximum  tangential  wind  (middle,  dashed),  and  the  radius  of  maximum  counterpropogating 
(outer,  dash-dot). 
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Fig.  10.  (a)  Relative  vorticity  field,  (b)  (9-field,  (c)  R-field  with  vorticity  contours  overlayed  to  relate  structures,  (d)  S-field,  (e)  integrated 
2-field  Q ,  and  (f)  FTLF-field  for  initial  time  of  tQ-lOh  with  integration  time  T- 1  h.  The  white  contours  in  (d)-(f)  mark  the  radius  of 
maximum  (inner,  solid),  the  radius  of  maximum  tangential  wind  (middle,  dashed),  and  the  radius  of  maximum  counterpropogating 
(outer,  dash-dot). 
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Fig.  11.  (a)  Relative  vorticity  field,  (b)  (9-field,  (c)  i?-field  with  vorticity  contours  overlayed  to  relate  structures,  (d)  S-field,  (e)  integrated 
2-field  Q,  and  (f)  FTLF-field  for  initial  time  of  to=12h  with  integration  time  T- 1  h.  The  white  contours  in  (d)-(f)  mark  the  radius  of 
maximum  ^  (inner,  solid),  the  radius  of  maximum  tangential  wind  (middle,  dashed),  and  the  radius  of  maximum  counterpropogating 
(outer,  dash-dot). 
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Fig.  12.  (a)  Averaged  angular  velocity  co(r,t),  and  (b)  averaged 
radial  derivative  dco/dr  for  r=0-48  h. 
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Fig.  13.  Radial  FTLE- values  a(ro,to)  for  integration  times  (a) 
30  min  and  (b)  1  h. 


rates  which  quantify  transport  across  shear-lines  (level  con¬ 
tours  Sc  of  S ,  see  Sect.  3.2.3). 

6.1  Measured  mixing  rate 

Given  a  level-contour  Sc  of  the  S-field  and  an  integration 
time  T,  we  define  the  mixing  rate  Rm  (C,  to)  as  the  area  of  ini¬ 
tial  conditions  whose  trajectories  cross  Sc  during  [to,to+T]9 
divided  by  the  length  of  the  contour.  This  mixing  rate  is  com¬ 
puted  (“measured”)  by  seeding  a  grid  of  initial  conditions 
and  counting  trajectories  which  cross  Sc- 

6.2  Predicted  mixing  rate 

We  define  a  predicted  mixing  rate  through  the  R -field  as  fol¬ 
lows.  Let  rc  be  the  average  radius  along  Sc-  If  R(xoTo)>0 
and  xo=ro(cos0o,sin0o)*  with  ro<rc,  xq  is  inside  the  circle 


with  radius  rc,  and  the  trajectory  0jo(xo)  is  repelling.  Thus 
trajectories  seeded  on  the  ray  with  angle  So  and  radial  val¬ 
ues  slightly  above  ro  move  outwards,  towards  the  circle 
with  radius  rc-  This  suggests  to  define  a  boundary  point 
(rc~8(Oo),Oo)  through  the  condition 

eR(-ro’eo’to)8(6o)+ro  =  rc.  (14) 

Points  on  the  Oo~ray  above  this  boundary  point  and  below 
rc  can  be  expected  to  cross  the  rc  -circle.  If  ro>rc  and 
R(xo,^o)>0,  trajectories  on  the  ray  0q  with  radial  values 
slightly  below  ro  move  inwards,  towards  the  circle  with 
radius  rc  again,  which  leads  to  the  same  boundary  point 
(Eq.  14),  now  with  <$($o)<0.  If  R(xo,fa)<0,  the  trajectory 
is  attracting,  and  initial  conditions  on  the  opposite  side  of  the 
rc -circle  move  towards  this  circle,  provided  xq  is  sufficiently 
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Fig.  14.  Radial  S-values  S(ro,to)  for  integration  times  (a)  30 min 
and  (b)  1  hr. 


Fig.  15.  Radial  /^-values  R(ro,to)  for  integration  times  (a)  30 min 
and  (b)  1  h. 


close  to  that  circle.  The  corresponding  boundary  point  is  then 
defined  by  (rc+8(6o),Oo),  where  8  satisfies 

eR(rc,e0,to)s (0O)  +  rc=  ro.  (15) 

By  varying  Oo,  the  conditions  (Eq.  14)  for  R>0  and  Eq.  (15) 
for  R<  0  define  an  annulus  of  initial  conditions  around  the 
rc  -circle,  whose  area  we  use  to  define  the  predicted  mix¬ 
ing  rate  Rp(C,to).  This  mixing  rate  is  an  approximation  of 
the  measured  mixing  rate.  Of  course,  several  approximations 
and  simplifying  assumptions  are  involved  in  this  definition, 
but  the  results  obtained  make  sense  and  the  structure  of  Rp 
is  similar  to  the  structure  of  Rm . 

Color-coded  plots  of  the  mixing  rates  a(ro,to),  S(ro,to)9 
R(r0,to),  Rm(C,to ),  and  Rp(C,to )  are  displayed  in  Figs.  13, 
14,  15,  16,  and  17,  respectively,  for  integration  times 
T =30  min  and  T= 1  h.  The  rates  a  and  S  should  be  compared 
to  the  average  angular  velocity  co  and  its  radial  derivative 


dco/dr  shown  in  Fig.  12.  These  rates  are  measures  of  shear 
and  show  the  highest  shear  during  the  initial  6h,  with  the 
amount  of  shear  dissipating  as  hyperbolic  mixing  begins  to 
occur  during  the  polygonal  eyewall  stage. 

The  hyperbolic  mixing  is  captured  by  R (ro,  to),  Rm  (C,  to), 
and  Rp(C,to).  All  of  these  rates  show  that  high  hyperbolic 
mixing  begins  with  the  polygonal  eyewall  formation  at  6h 
and  continues  through  the  transition  to  a  high  vorticity  eye  at 
24  h.  In  particular,  the  measured  and  predicted  mixing  rates 
Rm  (C,  to)  and  Rp(C,  to)  are  very  similar  in  structure,  and  re¬ 
veal  strong  mixing  near  the  zero  S-contour  So  (jet). 

We  note  that  Hendricks  and  Schubert  (2009)  studied  mix¬ 
ing  for  the  same  barotropic  model  using  the  concept  of  effec¬ 
tive  diffusivity  (Shuckburgh  and  Haynes,  2003;  Nakamura, 
1996).  This  quantity  was  computed  at  the  initial  time  £o=0, 
and  with  an  integration  time  t  =  T  (flow  map  0q)  varying 
over  the  full  duration  of  the  model  run  of  48  h.  The  effective 
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Fig.  16.  Measured  mixing  rate  versus  Rm(C,to)  across  shear-lines  Fig.  17.  Predicted  mixing  rate  Rp(CJq)  across  shear-lines  Sc  for 
Sc  for  integration  times  (a)  30  min  and  (b)  1  h.  integration  times  (a)  30  min  and  (b)  1  h. 


diffusivity  yields  a  mixing  rate  that  depends  on  r  and  t.  Our 
mixing  rates  have  a  similar  dependence  on  to  as  the  mixing 
rate  of  Hendricks  and  Schubert  (2009)  had  on  t ,  due  to  our 
moving  frame  approach.  The  mixing  rate  of  Hendricks  and 
Schubert  (2009)  shows  two  radial  regions,  at  approximately 
r- 30  km  and  r=50  km,  where  high  mixing  occurs  during  t=6 
and  £=32  h.  The  similar  time  and  spatial  regions  associated 
with  the  same  high  mixing  in  our  study  occur  for  a  and  S 
during  to=6  and  £o-24  h,  when  the  band  of  low  FTLE- values 
(the  jet)  is  less  prominent  and  more  hyperbolic  mixing  oc¬ 
curs.  The  mixing  rates  Rm(C,to )  and  Rp(C,to )  show  very 
similar  times  of  high  mixing  as  the  mixing  rate  of  Hendricks 
and  Schubert  (2009),  and  the  S-contours  for  which  these  high 
rates  occur  are  located  approximately  at  the  same  radial  val¬ 
ues  noted  above. 


High  R -values  determine  sets  of  trajectories  that  show 
growth  orthogonal  to  the  Lagrangian  velocity,  and  result  in 
filamentation  that  enables  turbulent  diffusion  to  occur.  The 
R -field  gives  an  advective  measure  that  converges  on  very 
short  time  scales,  yet  still  yields  similar  mixing  rates  as 
the  effective  diffusivity  obtained  by  integration  over  the  full 
model  time. 

7  Conclusions 

We  have  characterized  Lagrangian  mixing  in  a  two- 
dimensional,  nondivergent  barotropic  model  for  hurricane¬ 
like  vortices  through  several  diagnostic  techniques.  For  this 
model,  the  field  of  finite-time  Lyapunov  exponents  provided 
a  measure  of  total  particle  separation,  but  it  did  not  sepa¬ 
rate  the  effects  of  hyperbolicity  and  shear,  and  did  not  show 
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distinct  coherent  structures.  In  order  to  separate  the  effects 
of  the  high  shearing  occurring  in  the  model  from  hyperbolic 
mixing,  the  trajectory  separation  was  decomposed  in  direc¬ 
tions  along  and  normal  to  the  Lagrangian  velocity.  This  de¬ 
composition  gave  rise  to  two  Lagrangian  fields,  the  R -field 
and  the  S-field,  which  quantified  the  relative  contributions  of 
hyperbolicity  and  shear  to  the  mixing  process,  respectively. 
In  this  approach,  shear-lines  and  shear- strengths  were  identi¬ 
fied  with  level-contours  and  level- values  of  the  S-field. 

The  R -field  showed  coherent  structures  which  impacted 
the  mixing  and  mesovortex  interaction,  even  through  high 
shear.  The  outer  ridges  of  the  R -field  were  also  involved 
in  protecting  the  mesovortices  from  environmental  flow,  and 
the  persistance  of  the  ridges  was  associated  with  the  domi¬ 
nance  of  particular  mesovortices  over  others.  In  contrast  to 
the  other  Lagrangian  methods  showed  that  the  R-field  was 
able  to  distinguish  not  only  regions  of  high  mixing,  but  also 
the  structures  that  were  involved  in  the  evolution  and  annihi¬ 
lation  of  mesovortices.  Future  work  will  be  devoted  to  study 
the  impact  of  coherent  structures  on  mesovortex  interaction 
during  tropical  cyclogenesis  and  its  role  for  the  evolution  of 
a  wave  critical  layer. 

The  impact  of  the  coherent  structures  on  the  mixing  pro¬ 
cess  was  quantified  in  terms  of  time-  and  space-dependent 
mixing  rates,  with  the  spatial  dependence  displayed  as  a 
function  of  shear- strength  as  well  as  a  function  of  the  (av¬ 
erage)  radius  of  shear-lines.  Overall,  the  moving-frame  ap¬ 
proach  used  in  this  paper  provides  time-dependent  mixing 
rates  that  isolate  mixing  events  occurring  in  particular  time 
windows. 

The  methods  used  here  led  to  mixing  profiles  of  similar 
structure  as  in  Hendricks  and  Schubert  (2009),  but  with  a 
moving  time  frame  and  with  fast  convergence  in  Lagrangian 
fields  integrated  over  short  times.  The  S-field  provided  a  nat¬ 
ural  choice  of  contours  for  varying  initial  time  that  allowed 
to  quantify  mixing  by  determining  transport  across  them. 

Future  work  will  address  the  impact  of  hyperbolicity  in  the 
presence  of  high  shear  in  a  three-dimensional  setting,  where 
planar  shear  and  movement  orthogonal  to  the  shear  are  sep¬ 
arated  using  suitable  extensions  of  the  fields  introduced  in 
this  paper.  These  techniques  will  be  applied  to  a  realistic, 
three-dimensional  hurricane  model  that  has  both  shearing 
and  hyperbolic  components  governing  the  mixing  processes. 
The  focus  of  this  work  will  be  rather  different  from  that  of  a 
recent  paper  by  Sapsis  and  Haller  (2009),  who  studied  mix¬ 
ing  in  a  three-dimensional  hurricane  model  by  computing  La¬ 
grangian  quantities  along  a  slow  manifold. 
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